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Abstract 

We develop a stream function approach for the horizontal Hele-Shaw, Saffman- 
Taylor finger. The model yields a nonlinear time-dependent differential equa- 
tion. The finger widths derived from the equation are 1 > A > in units 
of half the width of the Hele-Shaw cell, in accordance with observation. The 
equation contains the correct dispersion relation for the creation of the finger 
instability. In an accompanying paper the stationary solutions of the equation 
are found numerically. 
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1 Introduction 



The Hele-Shaw cell experiment of a less viscous fluid displacing a more viscous one, 
is the paradigm of fingering phenomena[l]. Saffman and Taylor [2] found fingers that 
arise in the Hele-Shaw cell when oil or air or water, penetrate into oil or glycerin. 
An initial instability develops into a finger or several competing fingers. They com- 
pared the finger profiles to an expression derived from the analytical properties of the 
complex fiuid potential and the stream function. The agreement with the data was 
reasonable for large capillary numbers. 

Since then, a large body of works has added to our knowledge of the fingering 
phenomenon in various branches of the dynamics of continuous media, such as den- 
dritic growth, directional solidification, diffusion- initiated aggregation, flame front 
propagation, electromigration, as well as fingering in porous media. 

The fiow in porous media prompted the initial research of fingering phenomena. 
The topic is of the utmost importance in problems of transport in saturated and 
unsaturated soils, groundwater pollution, etc. ^ 

Despite considerable efforts, the fingering phenomenon remains in many aspects 
uncharted territory and it defies intuition. [4] As Tanveer[4] points out, small effects, 
like local inhomogeneities, thin film effects, and surface tension make the theoretical 
description very difficult. 

The original theoretical formulation of Saffman and Taylor [2] was ill-posed mathematically. [5] 
This is reflected in the indetermination for the asymptotic size of the finger. Con- 
trarily to the measured profiles, that were found to be bounded from below, by a size 
of around one half of the cell width, the analytical expressions showed no such lower 
bound. 

McLean and Saffman [6] improved the theoretical approach of Saffman and Taylor [2] 
by including the effect of surface tension. Tanveer[4] has shown that the approach of 
McLean and Saffman[6], is equivalent to an expansion in a parameter related to the 
finger half-width, capillary number and aspect ratio. This technique fails at the tail of 
the finger and higher order terms are needed. The results of McLean and Saffman[6] 
predicted profiles that matched very well the front (nose) part of the finger. They 
found a limit of A > | to the finger half-width. More recent experiments with dif- 
ferent aspect ratios of the Hele-Shaw cell found the finger width limit is more likely 
around A > 0.45 of the cell width[7]. 

There appear to be, other branches of solutions (two at least) [8] , that compete with 
the branch found by McLean and Saffman[6]. These branches are unstable[9, 10, 11]. 
Pitts [12] took advantage of the observed dependence of the curvature of the finger on 
the angle, to obtain an analytical scale covariant expression for the finger shape that 
fitted measured values (by Saffman and Taylor and by Pitts himself) extremely well, 
especially for finger half-widths smaller than A 0.8. For wider fingers, the profile 
function was found to miss the measured finger by a small amount. Pitts presented 

^The website http://www. maths. ox.ac.uk^howison/Hele-Shaw /helearticles.bib , cited in ref.[3] 
carries an extensive (more than 600 papers), up to date list of references on the Hele-Shaw problem. 
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analytical expressions for the finger profile including the effect of surface tension. He 
also derived an expression for the finger's asymptotic half-width containing both the 
capillary number, and a fitted parameter. This parameter takes into account the 
fluid fllms left behind by the passing flnger. Despite Pitts's success in reproducing 
the data, he did not provide a physical basis for the phenomenological assumptions. 

DeGregoria and Schwartz [13] used a boundary integral Cauchy technique to inves- 
tigate the production and propagation of fingers in the Hele-Shaw cell. Tryggvason 
and Aref[14] used a vortex in line method to determine the relationship between fin- 
ger half-width and fiow parameters. Both works, as well as the stationary calculation 
of McLean and Saffman[6] give similar results concerning the finger half- width de- 
pendence on capillary number. A marked improvement is found in the numerical 
calculations of Reinelt[15]. 

The present work offers an alternative path to the fingering problem. We discuss 
the validity of Darcian fiow in the Hele-Shaw cell. We find that the Hele-Shaw Darcian 
equation is valid in the direction of motion of the finger exclusively. Consequently, 
there is no harmonic condition on the hydrodynamic potentials that would prevent 
the existence of vortices at the finger boundary. 

From minimal assumptions and a convergent expansion of the stream function, we 
derive a nonlinear differential equation for the stationary finger. The equation is then 
generahzed to the time dependent case. We connect to the nonhnear hydrodynamic 
equation of Korteweg and deVries[16]. Numerical profiles for the stationary fingers 
and finger properties are displayed in the following paper. 

Section 2 discusses our departure from the conventional approach to the Hele- 
Shaw, Saffman- Taylor theory. Section 3 develops the stream function at the interface. 
Time independent and time dependent nonlinear equations for the finger, are obtained 
in section 4. Section 5 states our conclusions. 



2 Fluid Equations 

The dynamical creeping fiow equation for the Hele-Shaw cell is [1, 4]. 

where v is the viscous fiuid velocity in the plane of the cell, b is the thickness of 
the cell, fj, the viscosity, and p, the pressure in the fiuid. Eq.(l) is obtained by 
averaging the Navier-Stokes equations over the smallest dimension of the Hele-Shaw 
cell. However, the true and real problem is inherently three-dimensional [14], especially 
because variable width films of fiuid are left behind the advancing finger. The finger 
thickness is not constant comparing tail and nose. It is thinner at the nose and 
thicker at the tail. Reinelt[15] dealt with these films of fiuid by means of a three (or 
four) domain splitting of the Hele-Shaw cell. This method improved the agreement 
between theory and experiment. 
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Eq.(l) imitates Darcy's law for the flow in porous media[2]. Both displacing and 
displaced fluids are supposed to obey a Darcian law. Eq.(l) implies a vanishing 
vorticity and consequently a harmonic equation for either the velocity potential or 
the stream function. As noted previously in the literature [14] , the neglect of vorticity 
is unjustifled at the front. 

We here summarize the derivation of eq.(l) and conclude that it is mainly valid 
for the longitudinal direction of the cell. 

The Navier-Stokes equation in a horizontal cell for a Newtonian fluid[18] without 
gravity is given by 



p— = -Vp + piV'n (2) 

where = + u • V(ii) denotes the material derivative, p is the pressure and, jj, 

the viscosity coefficient. For the stationary case, there is no explicit time dependence. 
Also, higher order quadratic terms in the velocity are neglected, in the creeping flow 
Stokes approximation. The left hand side of eq.(2) then equals zero. 
Eq.(2) becomes 



Vp = yuV^u (3) 

The Hele-Shaw cell is usually a rectangular paralellepiped with a long dimension 
in the direction of the flow, the much smaller width, the y axis, and the 

smallest dimension by far, the thickness of the cell, the z axis. The aspect ratio 
of the cell, the ratio of the width along the y axis as compared to the thickness 
along the z axis is large. The cell is designed to be very thin. The most important 
contribution to the Laplacian in eq.(3) is the second derivative with respect to this 

small dimension, the z axis, and to a lesser undetermined extent, the — — r- term. The 

Poisseuille parabolic proflle flow with separable dependence of the velocities of the 
form (Eq.(4.8.19) in [18]) 



u = -—Vp{x,y)z{b-z) (4) 

solves eq.(3) with complete neglect of the derivatives with respect to x and y. Eq.(l) 
is then obtained by averaging eq. (4) with respect to the transverse dimension z. This 
brings in a factor of 6 in the denominator and a factor of 12 in the final equation (1). 

For eq.(4) to be valid, the pressure gradient has to extend to infinity[18]. This con- 
dition holds for the fiow x direction. However, it is not appropriate for the transverse 
direction. In the y direction there appear strictures comparable to the cell thickness 
between the finger and the lateral walls. Pressure gradients in the y direction can 
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be quite important. Even inertia effects may become relevant for the y component 
of the fluid equation of motion. The averaging procedure, that lead to eq.(l) is also 
known as the lubrication approximation. [18] In the lubrication approximation, even 
in the context of thin films[25], it is used for the direction of the flow only. The x 
axis is ihmited, whereas the y axis has a deflnite scale. That scale should be present 
in a correct treatment as a separate length scale, besides the thickness of the cell h. 
However, and due to the averaging procedure, it docs not. These arguments suggest 
that the transverse component of eq.(l) is not correct. 

In the following, we will use eq.(l) for the direction of motion of the flnger only. 
The flow in the transverse direction will be constrained by means of boundary con- 
ditions. 

The Darcian creeping flow equation we adhere to is 



. . d p . , 

where u is the velocity in the x direction, and p is the pressure averaged over the z 

coordinate. The equation for the velocity in the y direction will be left in the Navier- 

Stokes form (3). Contrarily to eq.(l), Eq.(5), does allow vorticity on the interface. 

Both eqs.(l,5) are not free of inconsistencies. [13] At the lateral solid edge of the 
w 

cell, y = ±— , with w, the width of the cell, we should demand impenetrability of the 
fluid and no-slip, namely 



IT = » 

OX 

Integrating the flrst of eq.(6) form the entrance of the channel to its exit, we 
flnd that the pressure is identical at both ends on the lateral side of the cell. The 
contradiction arises form the no-slip boundary condition applied to the creeping flow 
equation (5). This conundrum is usually resolved by ignoring the no-slip condition 
without further justiflcation.[13] Although, the use of no-slip boundary condition has 
been contested repeatedly, because of the non-integrable stress[19, 20, 21] it generates, 
^ and recent experimental work on the subject[23, 24] supports this view, it is not 
consistent to apply it for the the z direction and not for the y direction. Both directions 
are limited by the solid boundary of the cell. In the next section we propose a method 
for the partial restitution of the no-slip condition. 

Eq.(5) is supplemented by boundary conditions. 

At the interface, we have[15, 26, 4] 

possible theoretical resolution of the no-slip problem for a moving contact line may be found 
in ref.[22] 
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2 T T 
Ap = — C0S7 + — (7) 

where, T is the surface tension parameter of the Young-Laplace formula, R, the radius 
of curvature in the plane of the cell and 7, the contact angle between the finger and 
the cell wall in the transverse b direction. The first term in eq.(7) is usually dropped, 
because the flow depends on the gradient of the pressure. The term should be kept for 
three dimensional treatments that consider the variation of contact angle 7 and, the 
change of finger thickness as one proceeds from nose to tail. Pitts [12] circumvented 
the problem by introducing an effective parameter to account for the difference in 
pressure inside the fluid between tail and nose. The narrower the finger as compared 
to the channel width, the stronger the effect of films of fluid left behind [15]. In our 
minimal model we drop the first term in eq.(7). Therefore, our results are better 
suited for wide fingers and small capillary numbers. 
Equation (7) reduces then to 

Ap = ^ (8) 

with T an effective surface tension parameter. 

To complete the set of equations we assume the fluids to be immiscible. The 
velocity of the fluid at the interface has to be tangential to the boundary. As discussed 
in the next section, the immiscibility condition implies that the boundary is a stream- 
line. 

Finally, for an incompressible fluid of constant density, the continuity equation 
reads 

V • V = (9) 

The stream function ^, determines the fluid flow through v = Vx^. The 

— * 

stationary equation of continuity of eq.(9) is therefore obeyed by construction V • 
Vx* = 0. 



3 The stream function 

The interface between the displaced and displacing fluids in the Hele-Shaw cell carries 
vorticity[14]. Consider flgure 1, where a typical schematical proflle of a flnger in a 
Hele-Shaw cell is depicted. The tongue of less viscous fluid (a gas for instance) 
penetrates a more viscous fluid to its right. The picture shows a cross section or 
lateral view of a horizontal Hele-Shaw cell. The z coordinate is along the thin vertical 
direction. We have drawn a few vortices to guide the eye. 

Recent experimental works [27, 28], have shown that finger fiow involves vorticity. 
Gravity driven fingers advance by a kind of rolling motion, as suggested long ago by 
Yarnold[29], and West [30] in the context of capillary rise. 
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X 

Figure 1: Vortices on a transverse slice of a Gnger interface. 
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The all important property of the fluid flow, is not the irrotationahty, broken at the 

interface and elsewhere, but the solenoidality originating from the incompressibility 
condition. The stream function, is better suited for the theoretical analysis of the 
problem. 

We here determine the stream function at the interface by means of a convergent 
perurbation expansion. The stream function together with the Young-Laplace equa- 
tion yield nonlinear differential equations for stationary and time dependent fingers. 

X y 

We first renormalize the axes to x ^ x — — y ^ y — — 7- with w being the 

w/2 w/2 

width of the Hele-Shaw cell. As a consequence of this rescaling, the finger transverse 
dimension becomes less than one, and serves as an expansion parameter. 

Two conditions constrain the stream field: The impenetrability of the cell wall 
and the immiscibility of the fiuids. 

The corresponding stream function boundary conditions read 



ox 



(10) 



at y = ±1, and 



977 9^977 

dt dy dx dx 

ai y = f], where f]{x,t) denotes the interface curve. Eq.(ll) may be easily derived by 

drj di] 

writing the differential oi ri as dti = — dt + — — dx and using the definitions 

at ox 

_ dx 9* 
^ dt dy 

dr] 

For the stationary finger, in its rest frame we have 'q^ ~ ^- With this substitution, 
eq.(ll) becomes 

d"^ dri (9^ 

= aty = 1] (13 
oy ox ox 

The differential of ^ rf^ = — — dy + — — dx, together with eq.(13) determine the 

dy ox 

unique solution for the stream function in the stationary case at the interface to be 

^stationary = COnstaut (14) 
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Without loss of generality, this constant may be taken to be equal to zero, because 
the stream function enters the calculations only through its derivatives. 
An alternative expression of eq.(13) is 



a{v) = tan ^{r]') 



(15) 



where a is the angle between the fluid velocity and the axis of propagation of the 
finger (x axis) and rj' is the derivative to the finger (y axis) with respect to the x axis. 
For the time-dependent case, the solution of eq.(ll) reads 



Eq.(16) solves eq.(ll), with given by the appropriate solution of eq.(13). 
At the front, we can rewrite the eq.(16) in the form 



where rj{x,t) is understood as a function of x and t only. Eq.(17) shows that the value 
of \E' at the interface changes with time. This equation is needed for the development 
of the time dependent finger equation of motion. 

We now proceed to derive the stream function at the interface for the stationary 
case corresponding to a fully developed finger travehng at constant speed, to be 
generalized later to the nonstationary situation. 

Consider the stream function in the finger rest frame. The finger is traveling in 
the positive x direction at constant speed U. At long distances ahead of the finger 
the fiuid is assumed to fiow with a constant velocity V. As discussed in section 2, 
the no-slip condition cannot be imposed in the Hele-Shaw cell and consequently the 
velocity V is the same all over the cell width including at the lateral boundary. This 
is a reasonable assumption provided the viscosity of the fiuid is small. However, for 
large capillary numbers, or large viscosity, this procedure is unreasonable. Boundary 
layers next to surfaces are patent in common phenomena even for moderately viscous 
fiuids. 

We follow a more conservative approach and demand the stream line next to the 
lateral edge to carry a velocity V that is not the asymptotic velocity of the fiuid ahead 
of the finger. It depends on the fiow properties. For capillary numbers tending to 
infinity it has to be equal to zero. Large capillary numbers may be implemented by 
using very viscous fiuids. For such fiuids, the no-slip condition is a must. Therefore 
in the limit of infinite capillary number V—0. In the opposite limit, i.e. no viscosity 
at all (finger half-width going to the full cell width) it is equal to the finger velocity. 
Hence we can assert < V < U.^ The parameter V of our slip boundary condition 
will be determined selfconsistently by the dynamics, it is not a free parameter. 

^Reinelt's equations[15] carry an extra term at the lateral edge that transforms the velocity y to 
a variable. 




(16) 




(17) 
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Expanding the stream function to 0{y^) with a separable ansatz, antisymmetric 
in y for symmetric fingers, and in the finger rest frame, wc find (we drop from now 
on the tilde on x,y, i.e. we work with rescaled coordinates) 



y,t)^y {V{x, t) - U{x, t) + A{x, t) + B{x, t) + C(x, t) + ...) (18) 
with A, B, C unknown functions. 

Eq.(lO), as well as the constancy of the fluid velocities, for the stationary finger 
in its rest frame imply 



a* 

= ^= V-U ,at y = ±1 (19) 

Inserting eq.(19) in eq.(18), we are able to fix the functions B and C in terms of 
A. Exphcitly 

= A + 3S+5C 

= A' + S' + C" (20) 
primes denoting derivatives with respect to x. 

Recalling that constants are irrelevant for the stream function, the unique solution 
to this order becomes 



B^-2A 

C ^ A (21) 
The stream function of eq.(18) with the conditions of eqs. (20,21) reads 

*(a;, y) = {y-U)y + A{x) y {I - y'f + ... (22) 



with 



^ = aty^T] (23) 

In order to assess the convergence of the separable expansion, we consider the 
the stream function to 0{y^^). Such a high order is needed for the fulfillment of the 
boundary conditions. Proceeding in the same manner that lead to eq.(22), we find 



^i{x, y,t) = y [V{x, t) - U{x, t) + {A{x, t) + y' D{x, t)) (1 - y'f] (24) 



10 



with D{x, t) an unknown function. Only the x component of the velocity is needed 
in eq.(5). Including the higher order terms of eq.(24) this velocity becomes 



V, = V{x, t) - U{x, t) + {l- y'') {A{x, (1 - 5 y^) + D{x, t) (7 - 11 y'')) (25) 



However, for the stationary finger, we demanded = at the boundary (23), then 
D{x,t) = at the interface. Therefore, for the velocity at the interface in the longi- 
tudinal direction eq.(22) is exact. The equivalence between eqs. (22,24) breaks down 
for locations off the interface. Nevertheless, even inside the fluid, the y^ dependence 
with y < 1 makes the contribution of the higher order term unlikely. These terms 
correspond to higher order fluid velocity gradients. These gradients become relevant 
when there exist large stresses inside the fluid. However, eq.(5) was derived under 
the assumption of Poisseuille creeping flow, meaning an orderly stationary almost 
constant very slow flow with neghgible internal stresses in the plane of the cell. The 
large stress situation is beyond the scope of the problem at hand and higher order 
terms should be ignored even inside the fluid. 

In the next section we use eq.(5) to derive a differential equation for the stationary 
finger and another one for the time dependent situation. 

4 Nonlinear differential equations for the finger 

Inserting eq.(23) in eq.(22), A{x) is determined at the front to be 



Using eqs. (12,22,26), the velocity of the fluid at the static finger, the only place 
where we can determine A in closed form, becomes 



A{x) 



U-V 



at y = r){x) 



(26) 




'^x,static 



{V -U)+ A{x) (1 - rf) (1-5 rj^) 



"^y, static 



A{x) 




(27) 




with 77, a function of x. 

The rescaled curvature in 



eq.(8) is 



R (l+77'2)l 



(28) 
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Where primes denote derivatives with respect to x. 

The sign of the curvature is the appropriate one. This can be seen by using eq.(5). 
The left hand side of the equation is a the positive velocity of the finger, chosen here 
to move from left to right, therefore the curvature has to decrease along x. In the 
upper part of the finger the curvature is negative and becomes more so as we proceed 
along X. Therefore the right hand side of eq.(5) is positive with the positive sign in 
eq.(28). 

Eq.(28) is the standard expression for the curvature, and, it is easily derivable from 
the definition of the arclength and the corresponding angle for differential increments. 
Equations (5,8,28) after rescaling imply 



Tb^ d 

Vx = 



3 w^jj, dx 



(1 +77'2)f 



(29) 



The stationary finger equation is obtained now by evaluating the derivative in 
eq.(29), with 



Vx,static ^{V-U) (30) 



4 rf 

Transforming back to the rest frame of the cell, with the finger in motion, = Vx,static'\' 
U . we find 



V... 3 ^^^^ + (l + ryj^ 

(31) 



with 



W{r)) 



A erf + b rf 



\ — fj'i 

1 _ l2^lU 'Up- 

B ~ f ly' 
V 



e — 



U 



(32) 



and the suffix in eq.(31) indicating differentiation with respect to x. 

An alternative expression in terms of the arclength measured from the tail of the 
finger reads 

= n +!^^ + (l-r,')^ 
u ^^^^ + 1 _ ^2 ^ v-L '^Is) 4^ 

(33) 
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Customarily e of eq.(32) is equated to A, the finger half-width. In section (3) we 

discussed the no-slip condition in the Hele-Shaw cell and concluded that this proce- 

V 

dure is inconsistent. In the present work we leave e = — as a free parameter to be 

determined self-consistently from the solutions to the nonlinear equations. Asymp- 

5 — 1 

totically far back at the tail we must have W{r]) — 0, or, e = — —tt, — ■ A determines 



e. 



4 A2 



Eq.(31) implies a lower bound on the finger thickness. In the limit of infinite 
capillary number e — > and at the tail of the finger W{rj) = yields 



A> ^ 
a/5 



(34) 



This straightforward and simple result coincides with experimental observation 
and improves the lower bounds found in the literature by much more laborious means. 

The potential W guides the propagation of the finger. Eq.(31) is of third order in 
the spatial derivatives. Third order differential equations for interface propagation are 
well known in the literature. Some notorious examples are: Landau and Levich[31], 
Bretherton[32], and Park and Homsy[26]. 

The third order equation (31) can be transformed to a more manageable second 
order one in terms of the angle tangent to the curve, with 77 obtained by integration. 
For ds starting at the tail, where 9 ~ n, ds cos{9) — — dx. 

Equation (33) now reads 



d'O W{rj) 



77 = X — j ds sin{9) 



(35) 



With W{r]) defined in eq.(32). 

We now proceed to derive the time dependent equation. Using eq.(17), and re- 
calling the definition of in terms of the stream function of eq.(12), we now have 



-\- U -\- Vx,static 

Vx 



(36) 



w 

Rescaling as before the length coordinates 77, and x by — , and the time by 

t — —, and inserting eq.(36) in eq.(29), we obtain the time-dependent differen- 
tial equation 



= 



Vx 



1+77, 



4 B 



W(r)) - ^ 

Vx 



(37) 
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In terms of the arclength eq.(37) becomes 



Vss Vs 



1-77: 



4 B 



Vs 



(38) 



While for the angle tangent to the front, the time-dependent equation (37) reads 



d'9 1 



wi,) + cosie) -/- 



rj = X- Jds sin{e) (39) 

The time- dependent equation (37) gives the correct dispersion relation for the 
perturbation of a planar front as we now show. 

The dispersion relation we seek is[34] a = \a\ (1 — Xa^) with a the instability 
parameter, a related to the wavenumber of the perturbation and A a parameter pro- 
portional to the surface tension. The velocity, that is proportional to a, destabilizes 
and the surface tension stabilizes the perturbation. 

Consider eq.(37) at y ~ for a perturbation protruding from a moving front 



= ,..._3^ + (l + ,a»J^(l-^) (40) 



Eq.(40) can be rewritten as 
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4^ B—U^ +x - [dx^ (41) 



Following Chuoke et al.[34] we consider the perturbation 

x(y,t) = (3{a)e'^ 

= ay + at (42) 

P is the amplitude of the perturbation, a function of a[34]. The integral in eq.(41) 
may be rewritten as 

/ = y dx— = j Xy Xt dy at y = rj (43) 
Evaluating the integral (43) with eq.(42) we find 
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/ ^ ^ e''^' {e^'^y -1) 

^ a I3^a y (44) 
The curvature may be calculated now by using 



yy 



-Xyy 
,2 



= -a^ X (45) 

Using eq.(45), the integral of eq.(44), and, x of eq.(42), at i = to 0{y^), as we 
are assuming a perturbative expansion, x of eq.(41), becomes 



Xo ~ (-4 S q;2/3 + /3) (1 + ay) - a y (46) 

The constant piece of xo i^i eq.(46), —4 B a^j3 + (3, is irrelevant. This is the 
reason we do not need to specify the lower bound of the integral in eq.(43). We took 
it as y = 0. Any other choice will merely change the unimportant constant in x- 

To satisfy eq.(41), to the lowest order in y, eq.(46) has to obey the algebraic 
condition 



-4 B a^p + pa- ap^a = (47) 

The third term in eq.(47) is quadratic in p. Independence from the initial per- 
turbation amplitude, requires /9 to be a function of a, as assumed in eq.(42). We 
can now determine this functional dependence to be Pa = k, with k, a constant. 
P is a positive number, therefore, k has to be positive for positive a, or negative for 
negative a. 

Equivalently, the proportionality above can be written as 



P-ri- (48) 
\a\ 



Inserting eq.(48) into eq.(47) we find 



-4 B a^p + pa- apa^. = (49) 



\a\ 



Finally, redefining |q;| ^ 777, X — 4: B k^ we find the expected dispersion relation 

\k\ 
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cr = |q;| (1 — a a^) 



(50) 



the desired dispersion relation. 

The time-dependent equation was developed on the basis of an antisymmetric 
stream function, eq.(22), the solutions, may, nevertheless, lack any predetermined 
symmetry. Although the equations are symmetric under a certain transformation, the 
solutions can break the symmetry. Boundary conditions and initial conditions are key 
factors in determining the symmetries of the solutions. Such is the case of an initial 
asymmetric perturbation to the front. On the experimental side, a symmetric finger 
develops only asymptotically at long times. The transient behavior is asymmetric. 

5 Conclusions 

The contribution of the present work consists in the development of nonhnear differ- 
ential equations for the Hele-Shaw Saffman-Taylor fingering problem. The approach 
is based upon the treatment of the stream function as a separable potential without 
any free adjustable parameters. 

For the stationary case, the finger nonlinear equations are eqs.(31), and (33) the 
latter expressed in terms of the arclength, while eq.(35) is the corresponding integrod- 
ifferential equation, of second order, for the angle as a function of arclength. These 
are nonlinear and nonlocal equations. 

The time-dependent equation is of the third order in space, and first order in time. 
Eqs. (37,38), corresponds to r] in terms of x and the arclength respectively, whereas 
eq.(39) is its integrodifferential equation for the angle as a function of the arclength 
also. 

The equations found here, resemble the Korteweg-deVries (KdV) equation, that 
possesses solitary wave solutions[16]. The KdV equation describes the nonlinear prop- 
agation of shallow water waves in a channel. 

The KdV equation reads 

fh? 2Th\ 3 Co , r|t^ ^ 
^0 r^^^^ + TTT'^'^^ + ^^(^0 H ) = (51) 

where h, is the depth of the channel, cq = y/ g h is the speed of propagation of linear 
waves in the channel, g is the acceleration of gravity and rj denotes the soliton profile 
propagating in the channel along the x axis with velocity U = — — , in a fiuid 

Vx 

with density p and surface tension T. This is a third order equation nonlinear in 
Tj. The specific nonlinear term, and, the potential found in eq.(31), differ from the 
quadratic nonlinear term of the KdV equation (51). The nonlinear terms in eq.(31), 
and eq.(51) arise from the boundary conditions at the interface between the front 
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and the displaced fluid as well as the surface tension pressure drop. Eq.(31) however, 
carries also the information of the lateral flow constraint that is absent in eq.(51). 
Moreover, eq.(31) is based upon the creeping flow assumption of eq. (5), while eq.(51) 
does not assume such a restricted version of the Navier-Stokes equation. Therefore 
the nonlinear terms differ. 

Nevertheless, the equations are similar. Both are nonlinear of third order in the 
shape variable, and of first order in time. Both equations describe the propagation of 
stationary and stable shape fronts in a fluid medium. 

Links between soliton equations and the Hele-Shaw flnger appear in a work of 
Kadanoff [33] , in which varieties of Harry-Dym equations were found to be related to 
flnger development. The above equations belong to a broad class of Rux-like partial 
differential equations, widely used in the literature. In the context of Darcian flow in 
the Hele-Shaw cell, Goldstein et al.[35] studied instabilities and singularities by means 
of an equation of this type. The main difference between the equations we derived 
here, and the one of Goldstein et al. [35] , is, again, the appearance of a potential term 
W{r]) in eqs. (37,39). For the Hele-Shaw cell, the cell provides the potential for the 
propagation of the finger, while the pressure difference determined by the surface 
tension provides the dynamics. 
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